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We study the effect of the dissipation on the quantum phase slippage via the creation of "vortex 
ring" in charge density wave (CDW) systems. The dissipation is assumed to come from the interac- 
tion with the normal electron near and inside of the vortex core. We describe the CDW by extracted 
macroscopic degrees of freedom, that is, the CDW phase and the radius of the "vortex ring", assume 
the ohmic dissipation, and investigate the effect in the context of semiclassical approximation. The 
obtained results are discussed in comparison with experiments. It turns out that the effect of such 
a dissipation can be neglected in experiments. 



I. INTRODUCTION 

o . 

The sliding motion of a charge density wave (CDW) is one of the most interesting phenomena in low dimensional 
conductors. Under an electric field smaller than a certain threshold field ex, the CDW is pinned by impurities, 
lattice and so on in the sample, and the electric conductivity is small. Above ex, the CDW starts sliding motion 
and causes a sharp increase of the electric conductivity, the spontaneous current oscillation, which is known by the 
t-H \ name of the narrow band noise, and other various peculiar electric conduction phenomena jO . These phenomena are 

■ well described by the classical- mechanical treatment of the Fukuyama-Lee-Rice (FLR) model , which treats the 
CO | CDW as a deformable object under pinning potentials ||. 

Below £t, the creep motion of the CDW occurs with a help of thermal noise in the usual temperature range. Then 
^— ^ | the conductance is given as 

croce-^W (1) 

Os ! 

*» — , where Vb is the effective barrier height. In the lower temperature range, however, the quantum effect is expected to 

■ be important. In fact Zaitsev-Zotov observed in o-TaS3 a temperature independent nonlinear conduction below £t 
in the temperature range below 10 K, which is considered to be due to quantum mechanical tunneling of the CDW 

I i H . The experiment was carried out in the two-probe configuration. He also found that the following dependence of 
""O ' the conductance a on the electric field e, 
Ch ! 

acxe- (£o/£)2 , (2) 

J> . where Sq ~ 1.1 x 10 4 [V/m] in the range 3 x 10 3 [V/m] < e < 1 x 10 4 [V/m]. This form is different from the one due 
to the Zener-tunneling or the creation of kink-antikink pairs The above nonlinear conduction below et is 
observed only in thin samples. 

In the framework of the semiclassical theory, the CDW current due to the quantum tunneling is related to the decay 
rate of the metastable state. Then the conductance is expressed as 

a cx e^*. (3) 

Here Sb is the Euclidean action for the "bounce" solution, which is the least action solution, corresponding to the 
proper tunneling process pjJlOtl- Duan proposed the "vortex shell" in the 3 + 1 dimensional space as the "bounce" 
solution, and obtained the same electric field dependence as that in eq. (|J) ju| . The vortex shell represents the 
conversion process from the single-particle current to the CDW current or opposite by expanding the dislocation loop 
(or "vortex ring" in terms of the phase) in the perpendicular plane to the one-dimensional direction as shown in Fig.[lj. 
The motion of dislocation loop causes the phase slippage [jl2] 13 1. The magnitude of £o in cq. (^J) that Duan estimated 
is, however, about 10 3 times as large as that in the experiment. This is because he assumed the homogeneous CDW 
system, the quantum nucleation can occur wherever in the sample with the same probability, and the energy gain due 
to the creation of a vortex shell comes from the interaction of the CDW with the external electric field. Maki JlJ] 
noticed that the experiment was carried out in the two-probe configuration, and considered the strong pinning by the 
electric probes at both ends of the sample at first. He reset the distance between the probes to the FLR coherence 
length in the case of 3D weak pinning to get the agreement with the experiment. The reason of the replacement is, 
however, not clear. 
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The motion of CDW is expressed by the extracted macroscopic degrees of freedom in these work but the dissipation 
due to the interaction between them and the other degrees of freedom is neglected, which can affect the conductance 
due to quantum tunneling jl6|,[l7]], and is important in comparison with experiments. In this paper we investigate 
the effects of the ohmic dissipation on the quantum phase slippage according to Caldeira-Leggett theory [|l6Ul7| . We 
write down the "bounce" action in §[n| in the case without dissipation, and add to this the dissipation term coming 
from the interaction between the extracted macroscopic degree of freedom and the normal electrons in § III. We see 
that the evaluated dissipation term can modify the dependence of the conductance in the electric field. In §|IV 
discuss the result obtained in 



we 



III with that of experiments. It turns out that such a dissipation can be neglected in 



experiments. 



II. FORMALISM 



We start from the Euclidean action of the FLR model |2| ^,|18[, which is given as follows 
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where we choose z-axis as the one-dimensional direction. Here <j> is the phase of CDW, N$ = (irv-pl x l y ) , m* effective 
mass, e the applied electric field and / is the condensation density and / ~ 2 in the low temperature range. (v x , v y , vp) 
stand for the elements of Fermi velocity, and are given by (\/2t x l x /h, \/2t y l y /h, 2t z l z sin(i z fcp)/?l), where t x , t y , t z and 
l x , ly, l z are transfer integrals and lattice constants in each direction, respectively. V5 mp (r) is the impurity potential 
energy density. We introduce x' — {vF/v x )x,y' — (vp/vyjy, z' = z, and r' = v-pTy m/m* and obtain the isotropic 
form, eq. (||). 

In order to calculate the tunneling probability, we must specify the state of the system just before the tunneling. 
Duan considered the homogeneous CDW [|lj, which means there is no macroscopic deformation in the CDW. In this 
case, however, the estimated value of £o in eq. (|^) is too large. Maki considered the CDW pinned by strong pinning 
centers, which are, at first, assumed to be electric probes at both ends of the sample We put 5<fr = <fi — <fi m , where 
4> m is the phase of the metastable state, that is, the phase configuration before the tunneling. Then we have 



SS E - 5S® + 5S& + 8S eRc 
5S® =Ao f dV {d^f 



5SW - A I dV {2 % (60)] ■ 9 M m } 
4e 



SS e fi oid = A I d r J^ e5( t>, 



(6) 
(7) 

(8) 
(9) 



where Aq = Nq fhv x v y ^ m* / m/ Av-p , \x = (r', x', y', z'), and SSp represents the elastic and kinetic energy of the CDW 

in the case of vanishing cf> m , and SSp 1 ^ the contribution from the deformation of the CDW in the metastable state, 
and SS e fi e id is the contribution from electric field. 

In the case of the homogeneous CDW, SSp 1 * 1 is negligible. However, in the case of the CDW pinned by strong 
pinning centers, CDW deforms so large, which is expressed between two strong pinning centers as, 

€8 

0m = Z z ( z - L), 

TiVF 



(10) 
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where L is the distance between strong pinning centers, one of which locates at z = 0. In the above, we neglect the 
effect of impurities in the bulk, which is taken into account by slight modification of the amplitude of eq. (|l(]) jlB) . 
Then eq. (g) becomes 



SSI* 



-2A 



9(t> m 
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dz' 


z=z Q •> 
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(11) 



where we note the contribution from one strong pinning center. Here it is assumed that Sip varies sharply only in the 
vicinity of z — zq . Equation ( pi] ) gives much larger contribution to the Euclidian action than eq. (|)J) in a realistic case, 
so we can neglect the contribution which comes from SS e a e id om y which is taken into account in ref. 11. 

We can express the above action by employing the radius of the "vortex ring" at imaginary time r', -R(t'), as 
follows, 



6S E = 16tt 2 A 



dr'i?(r')yi + 
AneAo (L + 2Q 



dR (rQ 



e x 2tt / &t'R{t'Y 



(12) 



where is positive zero point of R (r'). In the case without dissipation, "vortex shell" is considered to have a sphere 
shape with a constant radius R. In this case SSe is given as 



SS K = AR 



ER 



where 



.4 



R 

l6n 2 A In—, 

Sz 



167T 2 eA (L + 2^) 

£j = — £. 



We optimize SSe about R, and get the bounce action, 



AA 3 
27 E 2 ' 



(13) 

(14) 
(15) 

(16) 



Then we have the conductance, 



a cx exp 



(eo/e)' 
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(17) 



Here superscript (e) means extremum solution. This form is the same as that observed in the experiment ||. In the 
case of the pure system and the two-probe method, L is the distance between the probes, and then the value of Eq is 
much smaller than that of the experimental result. 

The action we have evaluated just above does not include the degrees of freedom of the normal electrons though 
the motion of the vortex ring accompanies the conversion from the condensed electrons to the normal electrons or 
opposite. So we take into consideration the interaction between the normal electrons and the CDW, which causes the 
dissipation. The latter is described by the macroscopic degree of freedom, that is the phase and the Peierls gap in 
the present case. 



III. EFFECTS OF DISSIPATION 



Now, we investigate effects of the dissipation. The motion of the "vortex ring" accompanies the conversion of the 
normal electrons to the condensed ones or opposite. It is considered that there are many normal electrons in the 
vicinity and inside of the vortex core even in semiconducting CDW system, such as TaS3. It is difficult, however, 
to deal with the interaction between the CDW and the normal electrons microscopically |ll|. Here we consider the 
dissipation phenomenologically in accordance with the Cardeira-Leggett theory |l^,[l7j]. In this paper, we assume the 
ohmic dissipation which comes from the two types of interactions. One is caused by the interaction between the CDW 



3 



phase and the normal electrons in the vicinity of the vortex core. The other is due to the interaction between the 
CDW and the normal electrons in the vortex core. First we consider the former. The interaction is expressed by the 
CDW phase around the "vortex ring" , 8<j>, as follows, 



SS^ M = TL [ dx'dy'dz'dr'dr" [H ^ V '> T ' ] " S<t> { ?> V '> tl . (18) 



This action gives the frictional term r]5(f) in the equation of motion of cf>. Here 7 = y£x£y/£ z , r' — y x' 2 + y' 2 and 
5<fi (x 1 , y', z', t') is given by 

ir r' - R(t') 
5cj> (x\ y', z', t') = --- arctan 7 -±-L, (19) 

when jR (t) is the radius of the "vortex ring" at imaginary time r. We execute the partial integral about r' and t" 
and get 

SS^M = -2r/ 7 4 / dr'dr"dr'd^ — 



|z' 2 + [r'-R (r')] 2 } \z' 2 + [r' — R (r")] 2 } 



^^P.n^e r -/-,-«, ,20, 

where we introduce the step function to remove the self-interaction of the "vortex ring" . The part that satisfies 
t' + t" = gives large contribution to the integral because R (r) = R (— r) and the part that satisfies r' — t" = is 
excluded due to the step function (|r' — r"| — £ z ). We integrate eq.(pO|) with respect to r' and z' , but these integrals 
are complicated because their ranges depend on the sample size. Here, we assume one of the upper limits of the 
integral variables, r' or z' , is so large and can be set infinity. Then SS^ [5(f>] is evaluated as, 



/oo Ctrl l i 

dr'dr" [L A + CR(r')]^f r 



dr" 



xlnJ— j— Le(|T / -r"|-6). (21) 
When the upper limit of the integration about z' is large enough, L\ and C are given as follows, 

La~^, C~kA (22) 

27 Cz 

where Z s is the radius of the section of the sample. On the other hand, if that about r' is large enough, we have 

La~^, C^^ln^. (23) 
27 47 ^ z 

Next, we consider the latter interaction, that is, the interaction between the CDW and the normal electrons in the 
vortex core. In the vortex core, the gap vanishes and the Fermi surface is recovered. Since we cannot define the phase 
in the vortex core, it is necessary to describe the motion of the vortex ring in terms of the degree of freedom except 
for the phase. Here we can describe it in terms of its radius. Then the dissipation term which derives the friction 
term proportional to 2ttR (t') dR (r') /dr' in the classical limit plays an important role. We can raise the following 
term for instance, 

6& M = M! f dr > dT "2^ (r') mr>)-R^)f 
d 4tt J ( T / _ T ") 2 

where rj v is the viscosity coefficient per unit length. In fact this term gives the above friction term in the equation of 
motion of R(r). We execute the partial integral about r and r', and obtain 



/OO 
dr'dT"i?(r')' 
-OO 



dr' dr 1 ' 



xln |T , T \ q(\ t '-t"\-&. (25) 
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Thus, the total dissipation term of the Euclidean action is 
5S d [6<f>] 



f dT'dr" f 2r/ 7 4 L A + 2 V1 4 CR (r') 

J — OO V 



R(r') 
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x In ■ 



i6(|r' 



(26) 



In principle we can get the "bounce" solution by variating the whole action with respect to R{t). It is, however, 
difficult in general. So, we evaluate the above action in the weak and the strong damping regimes |l!||2(J. In the weak 
damping regime, we can approximate the shape of the "vortex shell" by a sphere. Then the dissipation term derives 
the term which is proportional to the surface of the "vortex shell" and that proportional to its volume. Therefore, A 
and E in eq. (jl3|) is expressed as 
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IT 
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Then the conductance is expressed as 

cr oc exp [— Sb/Fi] 

where 



e — ei 



(27) 
(28) 

(29) 
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The condition that the assumption of the weak damping is valid is expressed as, 

?77 4 La "C 16^4q and e » ei. 



(30) 
(31) 

(32) 



Comparisons of eq. (0) and eq. (29) indicates that the dissipation increases e in eq. (jl^) and decreases the applied 
field s. Thus, the tunneling probability is suppressed in two ways. It is to be noted that in the range s ~ Si the 
assumption of the weak damping and then eq. (2^) are not valid. 

In the case that rjj^Lj^ 3> 16 A and e -C £%, the system is in the strong damping region. In this region, the action 
becomes small when the part of R{t') parallel to r' axis becomes large [fL9| . Hence, we can assume the shape of the 
"vortex shell" by the cylinder. Then we have 



9 f\ „ \ Ta 8TT 2 eA (L + 2£ z ) , 
SS E = 16tt 2 A ( -R 2 d + 2R d T d ) In ^ —eR 2 d T d 



6. 



1 



2T d 



4t77 4 L a + 2777 4 Ci? d + -, v7 ^ d j i?3 In 
where 2Xd is the height of the cylinder, and R d is its radius. We optimize the action about T d and i?d, and get 
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ln-^- 



(33) 

(34) 
(35) 



Then the conductance is given as, 
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a ex exp 



(36) 



where 



_ 4vFA/M87rMo+47H 4 L A ) ik 
e(L + 2£ z ) ln 6 

4?i3tiF7 (2777C + ?7v) In -f— 



(e) 



£9 = 



e (L + 2Q 



(37) 



(38) 



Thus the logarithm of the conductance can be proportional to the reciprocal of the electric field cubed. As seen from 
eq. (|29|) and eq. (^7|), there exists the electric field e* crosses from the strong damping to the weak damping regimes. 
The magnitude of e* is the same order with that of E\. 



IV. CONCLUDING REMARKS 



In this section, we estimate the effects of the dissipation in experiments of the conductance. The value of each 
parameter is obtained in refs. 6, 14, 21. 

We first discuss the case that the dissipation comes from the interaction between the CDW phase <j) and the normal 
electrons in the vicinity of the vortex core. The coupling constant ij in this case is considered to be proportional to 
the normal electron density near the vortex and then to be expressed as rj — gp^, where g is a coupling constant per 
electron and p\ is the normal electron density near the vortex. We estimate the effects in NbSe3 at first, where Fermi 
surface survives even at low temperature and then is expected to have similar magnitude to that in the bulk, p\. 
The magnitude of gp\ can be evaluated from the bulk conductance under large electric field, where the dissipation 
for the sliding motion of CDW in this material mainly comes from the interaction between the CDW phase and the 
normal electron p^|. Then we obtain 

^~2xlO-V]xiA[m], (39) 
16^o 

^ 4ChVF 3xl(r 16 [V]x— ^— [V- 1 ]. (40) 



4eA (L + 2&)e ' e(L + 2£ z ) 

Since L A is irl s /2-f(~ 3 x l(T 6 [m]) or (L + 2£ Z _W47 2 (~ 2 x 10 _1 [m]), eq. (|9|) is much smaller than unity. On the 
other hand, C is the quantity of O (1), and eq. ( fiO|) is much smaller than unity, too. Therefore, we conclude the effect 
of such a dissipation is negligible in NbSe3. Next we try to estimate the same effects in o-TaS3. In this case we do 
not know the magnitude of 77. It is expected, however, that it is similar to or smaller than that in NbSe3. Hence the 
effect of such a dissipation is negligible even in this material. 

Next we consider the case that the dissipation comes from the interaction between the CDW and the normal 
electron in the vortex core. In this case the dissipation has quantitatively similar effect both in NbSe3 and TaS3. So 
we estimate it for NbSe3. The motion of the vortex accompanies the interior normal electrons, which enhance the 
vortex mass and cause the dissipation. The enhancement of mass is discussed in Appendix. The dissipation coefficient 
r? v is estimated as follows, 



Vv - ire pl c &V£xf y p res , (41) 

where p rcs is the residual resistivity, p™ is the electron density in the vortex core, respectively. From eq. (|3l|), E\ is 
given as follows, 

£1 - 2 x lfr 7 [V]/-L[m]. (42) 

This value is much smaller than £0 ~ 3 x 10 _1 [V]/L[m]. Then we conclude that the dissipation effect can be neglected 
also in this case both for NbSe3 and TaS3. 

In order to get the agreement of the magnitude of £0 with the experiment, we must set L ~ 3 x 10 -5 [m], which is 
much smaller than the sample length L s = 3.2 x 10~ 4 [m]. So we consider strong impurities in the sample may play 
the role of strong pinning centers. 
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Before the end of this section, we would like to discuss the reason why the nonlinear conduction below £t is 
observed only in thin samples |2q| . In order to observe the stationary current the CDW in the bulk as well as that 
near the strong pinning centers such as probes must move. Both motions are due to the quantum tunneling in the 
case discussed here. The quantum tunneling in the bulk is considered to be due to the kink-antikink pair creation [0. 
The nucleation rate T of kink-antikink pair in d dimension is given by 

r oc d-^/o". (43) 

The conductance of the system is proportional to the product of the tunneling rate near the strong pinning center, eg. 
eq. ( |l7j ) and that in the bulk, eq. (f|3|). When the sample is thin and its radius, l s , is smaller than the FLR coherence 
length in the perpendicular direction to the 1-dimensional axis, £ x ,y, we can regard the system as 1-dimensional in 
the bulk. It is to be noted that £ x ,y is much larger than R^ and then the system can be 3-dimensional for the the 
vortex ring. Then the tunneling in the bulk is easier than that near the strong pinning centers, and the conductance 
is determined by the latter, that is, the nucleation rate of the "vortex ring]\ However, when the sample becomes so 
thick that l s is larger than £ x ,y, the system is regarded as 3-dimensional [|24j. Since T is too small for d = 3, we can 
not observe tunneling current in this case. Hence we can conclude that the nonlinear conduction below £t is observed 
only in thin samples |25| . 

In this paper, we consider the effect of the dissipation on the CDW. We can discuss it on the spin density wave 
(SDW) in the same way. Moreover, it is concerned with the quantum phase slippage in superfluid 4 He and 
ultrathin superconducting wires p7|. 
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APPENDIX A: THE EFFECTIVE MASS OF THE VORTEX CORE 



We investigate the effective mass of the vortex core, m*. The effective mass coming from eq. (|T^) is obtained as 
follows. The term in the equation of motion with respect to R (r) from eq. (|l2|) is given by, 



SR{t) 



= I 8tt 2 A) In — 



167r 2 eA Le 



-R(t) 



8 2 R{t) 
8t 2 



R(t). 



(Al) 



We assume that \8R(t) /8t\ <C 1, which shows that the velocity of vortex core is much smaller than the phason 
velocity. We use this condition, and expand eq. (Al). Then we get the effective mass per unit length m* from the 
coefficient of d 2 R(r) jdr 2 as follows, 



47T An In ■ 



r (e)' 
' m 

IT' 



(A2) 



Next we take into account the contribution from the normal electrons in the vortex core to m*. We add the following 
kinetic term to the Euclidian action, 



SS VC ~ 7 vp 



dr'7T 2 pl%^R(r') 



„(dR{T') 



V dr' 



(A3) 



The effective mass coming from the above term is given as, 



V m* 



(A4) 
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The comparison of m* and m* c shows that m* c is always negligible. Then we can take into account only the dissipation 
term and neglect others which come form the motion of the normal electrons in the vortex core. 
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FIG. 1. The phase configuration in sample. The top of CDW, where Q-r + <j> = 2nn, is shown by solid lines in the plane 
perpendicular to the one-dimensional direction. The dots show the vortex cores. Moving vortex core, that is, expanding "vortex 
ring" , converts the condensed electrons to the normal electrons or opposite. 








FIG. 2. The vortex shell in 3 + 1 dimensional space. 
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